Incompleteness of Spin and Orbital Angular Momentum Separation for Light Beams 
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The angular momentum of a circularly polarized beam of light in vacuum is calculated and 
compared with the value of the same observable for a wave packet of light. While the packet has a 
finite transverse and longitudinal extent, the span of the beam along the direction of propagation is, 
in principle, infinite. This simple but subtle fact renders incomplete the textbook calculation of the 
separation in a spin part and in an orbital part of the angular momentum of a beam of light. We 
demonstrate that an extra spin-like part must be added in order to preserve the intuitive fact that 
the spin angular momentum of a circularly polarized beam of light containing A'^ photons cannot 
exceed the value ±7V times the reduced Planck constant. Finally, the impact of this extra term 
is quantified by means of two examples: a Laguerre-Gaussian and a Bessel beam, both circularly 
polarized. 



Introduction. — In his Treatise of Electricity and Mag- 
netism [1| , Maxwell pointed out that the electromagnetic 
field carries both energy and momentum, and that the 
momentum can have both linear and angular contribu- 
tions. Today it is common knowledge that the angu- 
lar momentum (AM) of light can be thought of having 
two contributions: a spin part associated to polariza- 
tion (firstly theoretically investigated by Poynting [2| and 
then experimentally demonstrated by Beth [3|), and an 
orbital part associated with the spatial distribution of the 
field, as first recognised by Darwin [4] . About sixty years 
later, thanks to the seminal work of Allen and Woerdman 
[51, the topic of AM of light experienced a new renais- 
sance, producing a considerable amount of literature on 
the subject, including theoretical discussions about the 
separation of AM in spin and orbital parts for optical 
beams [6l-[li|. the spin-Hall effect of light [l5l-[l8l|. the ge- 
ometric spin- Hall effect of light 1£, 2^ , and a huge vari- 
ety of experimental applications such as optical tweezers, 
spatial light modulators and vortex beams. A satisfac- 
tory review of the developm ent of t he field of optical AM 
of light is reported in Refs. 



21-2 



In this Letter we report on a novel separation of the to- 
tal AM of a light beam into three parts: the orbital part, 
the spin part and a new spin-like part which arises from 
an incomplete two-dimensional integration (over a plane 
orthogonal to the direction of propagation of the beam) 
of a three-divergence term. The latter part, also called 
"surface term" [24|, is usually neglected when the light 
field has a finite extent along any direction. This is not 
the case for a beam, which is unbounded along the axis 
of propagation. We refer to this additional term as the 
surface spin angular momentum (SSAM). We show that 
SSAM arises whenever the intensity of a beam is recorded 
by a planar detector which, in practice, performs a two 
dimensional integration over a cross section of the beam 
perpendicular to its direction of propagation. As this is 



the case occurring in the majority of experimental opti- 
cal setups, we believe that the impact of our new results 
may be significant. 

The structure of this Letter is as follows. We first cal- 
culate the total AM density j(r) of an arbitrary electro- 
magnetic wave field. Then, by means of some straight- 
forward manipulations, we decompose j(r) into the sum 
of three terms: the well-known orbital and spin AM 
parts and a third term (the SSAM) which amounts to 
a three-divergence. Next, we integrate j(r) upon a two- 
dimensional surface at z = const., mimicking the de- 
tection process, thus obtaining the total AM per unit 
length. This unravels the nature of the SSAM term and 
how it manifests in the actual detection of the AM of 
a light beam. Finally, we calculate the explicit form of 
the SSAM for the cases of Laguerre-Gaussian and Bessel 
beams. 

Theory of SSAM. — To begin with, by following van Enk 
and Nienhuis [8| we write the cycle-averaged linear and 
AM vector densities in vacuo of a monochromatic field of 
angular frequency wq, as 

p(r) = |Re (E* x B) , j(r) = |r x Re (E* x B) , (1) 

where the transverse electric field E and the magnetic 
field B are expressible via the gauge invariant transverse 
part A of the vector potential as 



E = iojoA 



and 



B^ Vx A, 



(2) 



respectively. The complex-valued time-independent 
fields A, E and B, are functions of the position vector 
r solely. The corresponding time-dependent real-valued 
fields are given by A.{r,t) — Re(Ae^"^"*), et cetera. The 
total linear and angular momenta of the electromagnetic 
field are given by the integrals of the corresponding den- 
sities over all space: 



p(r)d3 



J 



/J<- 



(3) 



By using Eqs. ^H^ and choosing, throughout this 
Letter, the physical units in such a way that eQUJo/2 = 1, 
we can express the AM density as 

j(r) ^ - RejiA* • (r x V) A + A* x A 

-5,[^;(rxA)]}, (4) 

where, according to [2J|, the first term represents the or- 
bital part of the AM and the second term is the intrinsic 
or spin AM of light. As it is also noticed in Ref. [24J, 
the third term is a three-divergence that vanishes when 
integrated over a volume whose surface boundary is so 
far from the field sources such that A(r) sa 0. Thus, if 
we substitute Eq. ^ in the rightmost expression for JT" 
in Eq. ^ , we obtain the textbook expression for the or- 
bital and spin AM separation. However, as pointed out 
by Haus 25[ , for a monochromatic beam of light propa- 
gating in the z direction the classical optics analogous of 
the helicity of a photon, namely the projection of the spin 
AM along the direction of the linear momentum, is given 
hy Sz — z ■ 3 where, henceforth, the over line denotes 
normalization according to the rule 



(5) 



U/ujo' 



and with J and U we have denoted the cycle averaged 
AM and energy per unit length: 



J= /jWd'i?, U^iJof 



A(r)f 



■'d'R. 



(6) 



In the expressions above d R = dxdy and the integration 
is over all the plane xy at z = const. This definition is also 
consistent with the physical picture of a beam of light as 
observed on a laboratory bench, where it travels across 
plates, lenses and other optical devices before eventually 
shining the surface of a planar detector, as opposed to a 
3D particle detector. 

As a consequence of the surface, as opposed to vol- 
ume integration, a beam-like field "feels" the contribution 
given by the z-part of the latter term in Eq. (j4]) , namely 
dz[Al{r x A)] because such term is not integrated when 
one calculates the AM per unit length, namely the in- 
tegral of the AM density j(r) over all the xy plane at 
z = const. In practice, since J dj[A*{r x A)]d^i? = 

dz / A*(r X A)d^i?, we can write the SSAM contribution 
in the following form: 



rsurf 



d_ 

dz 



/"Re|iA*(r X A)|d^i?. 



(7) 



From Eq. ([7]) it seems that the 3^^'^^ depends both on the 
potential vector A and the position vector r, thus reveal- 
ing a seemingly hybrid orbital/spin nature. However, in 
the remainder we will show results that do not support 
this interpretation. 



Now, in order to give the SSAM an explicit form, we 
use the angular spectrum representation for the vector 
potential using the so-called helicity basis [26| ho-(K) = 
(l/V2)[ei(K) + i(Te2(K)], where ei^2(K) are two unit 
vectors orthogonal to k = K -|- fc^z, where we have 
written the transverse part K of the wave vector k as 
K — k^K + kyY and k'^ + ky + kl — k^ with fco = wo/c. 
Here with a = -\-l and cr = — 1 we denote left and right 
circular polarization, respectively. The helicity basis vec- 
tors ei_2(K) together with k = k/fco itself, form a "local" 
Cartesian reference frame {ei(K), e2(K),k} attached to 
the wave vector k, whose expression can be obtained from 
the "global" reference frame attached to the central di- 
rection of propagation z of the beam and denoted with 
{x,y, z}, via a rotation by the angle — arccos(k • z) 
around the axis n = z x k/|z x k| [27|. Explicitly, the 
potential vector in the angular spectrum representation 
reads as: 



A(r) 



-y 



a„{K.)h.^{K)e'^'''e^-^d^K, (8) 



where d A' = dk^dky and kz 



— \/kQ k^ ky > 



for the case of forward homogeneous fields con- 
sidered here [26|. Substituting this expression into 
Eq. d?]) and using the relation Re{jA*(rxA)} = 
i [z{A* X A) - r • (A* x A)z] /2 we obtain after a 
straightforward calculation: 



^■"■' = 5E'/km'(I:-^)^'- 



(9) 



This explicit form for J'^^''^ is the first main result of this 
Letter: What does it represent? The transverse part of 
the term under the sign of integration is equal to —1/2 
times the standard AM density in Fourier space, namely 
K/fcg- Conversely, the longitudinal component is a com- 
pletely new term that can be written as k^/(2v'1 — k^), 
where k = \K.\/ko. This contribution is O(k^) and there- 
fore disappears in the standard paraxial limit where such 
terms are not retained [25|. In order to understand how 
this surface part affects the total AM per unit length we 
must add it to the standard spin AM density originating 
from the second term in Eq. ^ and equal to 






(10) 



By comparing Eqs. © and (|10p . we immediately recog- 
nize that 3^™^ has the structure of a spin angular mo- 
mentum, as anticipated in the discussion after Eq. ([7|). 
Finally, the total spin AM is simply given by the sum 



fspm _ Tspii 
"J tot — "* 



Tsurf . 



Trspm 
"^tot 



<T = ±1 



l«.(K)|^ 



k 



kz^ 



d^K. 



(11) 



In the paraxial limit the terms between parenthesis in 
the equation above reduce to 2z and we recover the well 
known paraxial result Jtot" = ^J^a'^ I |fl(T(K)pd^if. If 
we denote with 9 = arccos(k • z) the angle formed by 
the generic wave vector k with the axis z, then the lon- 
gitudinal contribution to Eq. (Ilip can be written as 
z(cos 9+1/ cos0)/2. Interestingly, this term has the same 
form of the intrinsic AM with respect to the z-axis of a 

beam when observed in a reference frame tilted by the an- 

I — 1 1 — L 
gle 9 with respect to the direction of propagation |23I.I28|. 

For the sake of definiteness, in the remainder we as- 
sume that the beam is circularly polarized along the 
direction f — (l/\/2)(x + iSy), where E — ±1 is the 
configuration-space helicity assigned in the (global) lab- 
oratory reference frame {x, y,z} and not to be confused 
with the momentum-space helicity a defined with re- 
spect to the (local) wave vector frame {ei(K), e2(K), k}. 
Therefore we can write the angular spectrum of the beam 
as ao-(K) — a(K)acr(K), where ao-(K) — f •hcr(K) is sim- 
ply the projection of the global polarization of the beam 
onto the local helicity basis. From this definition it fol- 
lows that J2a\'^'^\^ = (l + '«2)/2 and X]ct^I"<tP = '^zS, 
where we have defined k^x -|- Hyy + k^z = k = k/fco- 
Thus, according to Eqs. ^H^ we can write 



fsurf 



and 



U 

UJQ 






(12a) 
(12b) 



(13) 



Since 3^^'^^ c>c S, then irrespective of the form of a(K), 
jsurf (jj-ops to zero for a non-circularly polarized beam 
thus revealing the spin, as opposed to orbital, character 
of the SSAM. From Eqs. (|12p we can straightforwardly 
calculate the projection of the spin and surface AM onto 
the central direction of proagation z of the beam: 



jspin.^^ J. f \a{K)\^nld^K. 
From these expressions it follows that 

pspin ^ jsurf) . £ = I /" |a(K)|2(l + nj) d" K 
Wo 



(14a) 
(14b) 



(15) 



where Eq. P^ has been used. Now, by using the defini- 
tion Eq. ([5]), we obtain 



Equation (1161) shows that the role of the SSAM term is to 
ensure that the z-component of spin AM per unit length 
along the propagation axis z is exactly equal to S times 
the energy per unit length of the beam, the latter being 
completely arbitrary. This our second main result. 

Examples. — In the remainder of this Letter, we calcu- 
late explicitly Eq. p^ for Laguerre-Gaussian and Bessel 



beams to investigate whether orbital AM affects J^^urf or 
not. Consider first the case of a Laguerre-Gaussian beam 
of waist Wo, angular aperture 9o = 2/(koWo) and with 
azimuthal and radial mode indexes £ and p, respectively 
[29J. In this case, from Eqs. (|12m3p we obtain 



Tsurf — 



/? 



1-/3 



jspin = S^-^Z, (17) 



where /3 = (6*0/2)2 (2p -f \(.\ + 1). Although P depends on 
\£\ such a dependence is not a consequence of the or- 
bital AM of the beam (this would lead to a dependence 
from the signed quantity £) , but it comes from the radial 
dependence of the beam profile. This result should be 
compared with the ones presented in Ref. [6[ . Note that 
for a fundamental Gaussian beam ^, p — >• and (3 — >■ 6'§/4. 
Next, we consider the case of a ^th- 
order Bessel beam of the form £'(r, 0, z) = 
J^(rfco sini9o) exp (i^0) exp (izfco cosi?o), where -do is 
the angular aperture of the characteristic Bessel cone. 
A straightforward calculation shows that Eqs. ([T7|) are 
still valid providing that /? = sin^ t?o/2. Here /3 is inde- 
pendent from the orbital AM value £, thus supporting 
the claim that the SSAM has a spin nature. Note that 
for both Laguerre-Gaussian and Bessel paraxial beams 
with 6*0 = i9o < 1, one has P^ - 0(6*^), namely J™^ 
becomes negligible. 
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FIG. 1. SSAM for a circularly polarized fundamental Gaus- 
sian beam(green line) and Bessel beam (blue line). 



Th e absolute value of the "surface helicity" I]'^"''^ = 
I Jsurf .g| ^ is depicted in Fig. [T]for both the case of a funda- 
mental Gaussian beam (green solid line) and of a Bessel 
beam (blue solid line) for i?o = ^0 € [0,7r/2]. From this 



figure it appears clearly that in the paraxial domain S''"''^ 
approaches zero and its contribution can be neglected for 
^0 ^ 1 ■ This is in complete agreement with the fact that 
in the paraxial limit /c^ ~ fco, K ~ and from Eq. Q it 
follows that J'^'^'f — >. 0. This term, however, plays an im- 
portant role in the non paraxial case, where, depending 
on the beam divergence, it can even become comparable 
with the standard helicity value |E| = 1 [,25|. 

Conclusions. — By means of a careful analysis of the 
standard procedure leading to the separation in the spin 
and the orbital parts of the AM of a light beam, we 
have shown that a self-consistent derivation of this result 
should account for the fact that the AM per unit length, 
as opposed to the total AM, is the physical quantity actu- 
ally measured in most experimental optical setups. This 
practically means that thoroughly all the steps of such 
derivation 3D integrations over the whole space must be 
replaced with 2D integrations over the plane xy at z = 
const. As a result of this replacement we have found that 
a new surface spin-like part, the SSAM, must be added 
to the conventional spin and orbital terms. Further, we 
have demonstrated that the SSAM term guarantees that 
the z-component of spin AM per unit length along the 
propagation axis z equals the helicity S(= ±1) times the 
energy per unit length of the beam, and that this result 
is universal, namely valid for any beam either paraxial 
or non paraxial. Finally, by calculating the SSAM for 
the cases of Laguerre-Gaussian and Bessel beams, we 
have shown that this novel term may give a contribution 
comparable to the conventional spin angular momentum. 
We believe that our results may have a relevant concep- 
tual impact upon the very lively and timely research field 
about AM of light [Ulsil. 
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